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ABSTRACT: We argue that the field theory that descibes randomly branched poly- 
mers is not generally conformally invariant in two dimensions at its critical point. In 
particular, we show (i) that the most natural formulation of conformal invariance for ran- 
domly branched polymers leads to inconsistencies; (ii) that the free field theory obtained 
by setting the potential equal to zero in the branched polymer field theory is not even clas- 
sically conformally invariant; and (iii) that numerical enumerations of the exponent 9{a), 
defined by Tj^{a) ~ \^ N^^i'^)^'^ ^ where T]\j{a) is number of distinct configuratations of 
a branched polymer rooted near the apex of a cone with apex angel a, indicate that 9{a) 
is not linear in l/o; contrary to what conformal invariance leads one to expect. 
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I. Introduction 

It is widely believed that statistical mechanical systems at the critical point of a second 
order phase transition are conformally invariant on scales much larger than any micro- 
scopic distance^^^^"''''^".'^ In two dimensions, the conformal algebra is infinite dimensional 
and conformal invariance places strong constraints on the fixed point correlation func- 
tions. The analysis of these constraints, starting with the paper of Belavin, Polyakov, 
and Zamolodchikov^^^'*^^'^, has lead to a rather complete understanding of two dimen- 
sional critical theories^'^''^^. The conformal properties of most random geometrical sys- 
tems: linear polymers, branched polymers with fixed topology, theta polymers, percolation, 
dense polymers, etc. in two dimensions are well understood'^"^^'^'^**^^. Randomly branched 
polymers^^^^"'*'''^ are a notable exception. Numerical enumerations of randomly branched 
polymers in the wedge^^^^"' , and transfer matrix calculations of randomly branched poly- 
mers in the strip*^^''^*'^"^, yield results at variance with the naive predictions of conformal 
invariance. 

In this paper we argue that the critical point of the theory describing randomly 
branched polymers is not generally conformally invariant in two dimensions. The ar- 
guments are quite simple. We state them here in outline. 

(i) We first show that the most natural formulation of conformal invariance for the field 
theory that describes randomly branched polymers^"^^"'*^^'^"''**^ leads to inconsistencies. 
We proceed by first assuming the field theory in question to be conformally invariant, and 
then by studying the consequences of this assumption. 

In order to study the consequences of conformal invariance, we have first to say how the 
fields in the randomly branched polymer field theory transform under conformal transfor- 
mations. This leads us to "the most natural formulation" of conformal invariance for ran- 

More precisely, translational, rotational, and scale invariance, in systems with only short 
range interactions, is thought to imply conformal invariance. 
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domly branched polymers. The most natural formulation of conformal invariance follows 
immediately from the supposition that the most relevant fields in the randomly branched 
polymer problem are those which appear in the effective Hamiltonian that defines the field 
theory (see equation 1 below). This is a very mild requirement; the field theory in question 
describes the large distance physics of randomly branched polymers, so we expect that the 
most relevant operators in the branched polymer problem to be those of the field theory, 
i.e. those which appear in the effective Hamiltonian, and which are integrated over in the 
functional integrals. 

Given, as seems likely, that the most relevant operators in the randomly branched 
polymer problem are those appearing in the effective Hamiltonian, it follows that the most 
relevant field in the effective Hamiltonian, the field with lowest scaling dimension, is pri- 
mary in the sense of conformal field theory (that is to say, if the theory is conformally 
invariant, then the lowest dimension field must be primary). For if this field were not 
primary, it would, by the usual structure of conformal field theories, have to be a sec- 
ondary field, and therefore the descendent of a field with smaller scaling dimension, and 
by assumption there is no such field. 

Given that the lowest dimension field in the effective Hamiltonian is primary, it is 
staightforward to show, by conformally mapping the plane to the cone, that one obtains 
inconsistent predictions for the exponent d{a) defined by Tjq{a) ~ N~^^'^^'^^ , where 
T]\j(a) is number of distinct configurations of a branched polymer rooted near the apex of a 
cone with apex angle a. From this we can conclude that either the most relevant operators 
in the branched polymer theory do not appear in the effective Hamiltonian (which, we 
repeat, seems to us unlikely), or that the field theory is not conformally invariant. 

(ii) The field theory that describes randomly branched polymers has a Parisi-Sourlas 
supersymmetry i''^''***. To further investigate the question of the conformal invariance of 
randomly branched polymers, we consider the simplest theory that has the same kind of 
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Parisi-Sourlas supersymmetry. This simple theory is obtained by setting the potential 
equal to zero in the branched polymer effective Hamiltonian. We show that this free field 
theory is not even classically conformally invariant. The reason for this is that the free 
Parisi-Sourlas theory is essentially a higher derivative Gaussian theory, and it is known 
that, in higher derivative theories, scale invariance does not imply conformal invariance 
zznn r^Yie lack of conformal invariance in the free theory does not of itself imply that the 
interacting theory of interest is not conformally invariant; it does however make it seem 
less likely. 

(iii) Finally, we present numerical enumerations of of randomly branched polymers in 
the cone. If randomly branched polymers were conformally invariant, the exponent 9{a) 
in the cone would be simply related to the same exponent in the plane. More precisely, 
conformal invariance leads one to expect that 9{a) should be linear in ^. Our enumerations 
indicate that this expected linear relation does not hold. 

The plan of this paper is as follows. In section 2 we present the natural formulation 
of conformal invariance for randomly branched polymers, in section 3 we show that this 
formulation leads to inconsistencies, in section 4 we show that the free Parisi-Sourlas 
theory is not conformally invariant, in section 5 we present our results for the exponent 
9 (a) obtained by exact ennumerations of branched polymers in the cone, and in section 6 
we restate our conclusions. 

II. Natural Formulation of Conformal Invariance 

In this section we present the most natural formulation of conformal invariance for 
randomly branched polymers. The effective Hamiltonian of the field theory that describes 
the universal properties of randomly branched polymers^^^^"'*^^ can be written^'^''*^* 

Hb.P = J A[a;(-V20 + F'(0)) + ia;2 + v;(-V2 + F"(0))V.]. (1) 

where 1^(0) = ^rcfP' + and lj are commuting scalars; and t/S anti-commute. p 
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is invariant under the supersymmetry transformations S(f) = — ae^x^V^, Sco = 2ae^5^7/', 
Sijj = 0, and Sip = a{e^x^u — 2ef^d^^(j)) where a is an anti-commuting number and tf^ an 
arbitrary vector^""''^"'*. 

The supersymmetry imposes relations between the the correlation functions of the 
theory. For our purposes, the following idenities are useful 

(0(r)a;(O)) = (V'(r)V'(O)) = A^{cj>{r)m)- (2) 

From these equations it follows that the the scaling dimensions a;^, x^j-, x^, and x^ of the 
fields (f), uj, '4>i and t/S are related by 

^(f) — ^u) — 2 = x,fp — 1 = x^ — 1. (3) 

(j) is therefore the lowest dimension field in Hq p . 

As stated in the introduction, the most natural formulation of conformal invariance for 
randomly branched polymers follows from the supposition that the most relevant fields in 
the randomly branched polymer problem appear in the effective Hamiltonian eq. 1. If this 
is the case, and it would be quite unusual if it were not, then it follows, by the argument 
given in the introduction, that the lowest dimension field, 0, is primary. 

Given that is primary, it follows from equation (2) that 

u = kV'^cj) + u (4) 

where a) is a sum of primary operators, and the constant k = This is because in a 

conformal field theory the two point function of operators belonging to different conformal 
families is always zero; since (0a;) is not zero, u> must have a term that belongs to the 
conformal family of (f). Since the scaling dimension of a; is a;^ + 2, this term must be a 
descendent of at level two, and the only descendent of (f) at level two which is also a 
scalar is V^0. Likewise, u) must be the sum of primary operators, since if o) has a term 
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which is not primary, that term must be the descendent of another operator with scahng 
dimension where n is a positive integer. But again there is no such operator in 

the effective Hamiltonian and therefore, in accordance with our assumptions, none in the 
theory (other than 0; a) is a commuting variable and therefore cannot be a first generation 
descendent of or i/;). It is important to note that conformal invariance together with cq. 
2 imply that uj = k'V'^cf) + uj holds as an operator identity, and not just as an equality in 
the plane. 

The statement that is primary, together with eq. 4, constitute our formulation of 
conformal invariance for branched polymers. 

III. Consequences of Conformal Invariance 

We now show that this natural formulation of conformal invariance for randomly 
branched polymers leads to inconsistencies. 

The correlation function (0(r)a;(O)) is the generating function for the number wj\i{r) 

of randomly branched polymer configurations with N bonds containing the sites and 
„shapir,miller 



Here is a fugacity for the number of bonds in the polymer, Kc is the (non-universal) 
value of the fugacity at which the average polymer size in the fixed fugacity ensemble 
diverges, and e measures the deviations of K from Kc'. e = ^"^^ . Equation (5) holds 
for e small, and ~ means that both sides of the equation have the same leading singular 
behavior. In the scaling region we expect {(f){r) (f)(0)) = , where the correlation length 

^ ~ e~'^. From equation (2), {(j)(r)uj{0)) = ^il^+l • Using this scaling form and integrating 
both sides of equation (5) with respect to r and then taking the inverse Laplace transform 
with respect to e we have 




(5) 



t: 



-lunrooted 
N 




N j^i^{d-2x4,-2)-3 



(6) 
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where Tj^p'''^^^^^ is the number of unrooted branched polymers of size N in the plane, and 
Z^r '^-/v('') — iV^T^"^'"'"'*^'^. Setting d = 2 and using the approximate value of = ^Q^derrida 
and the exact value of ^ = 1 P"'"^'^^ we find = —l/u = —1.5. The field (j) therefore has 
negative scaling dimension. 

Using our proposal for conformal invariance, we can also calculate the exponent 9{a), 
defined by T/v(a) ~ N-^i<^)+'^^ where T/v(a;) is the the total number of randomly 
branched polymer configurations rooted near the apex of a cone with apex angle a.^ To 
calculate 9{a) we need to know (0a;) in the cone. The transformation 2; — > C = maps 
the plane onto the cone. Since, by assumption, (f) is primary, we have at criticality 

{<P{Z1)<P{Z2)) plane = l§(^l)n^(^2)r^(</>(Cl)</'(C2))cone (7) 



SO that. 



mi)m))cone = (-)^^^|Cl|(^-^^"lC2|(^-^)-^ \. (8) 



where we have normalized so that {4>{f)4'{^)) plane = • f^om eq. 4 it then follows 
that 

(0(ClV(C2))cone = (0(Cl) + ^))cone = (0(Cl)0(C2))cone (9) 

where we used the orthogonality of (j) and a) to obtain the second equality. Combining eqs. 
8 and 9 then gives 

The correlation function (0(Cl)<^(C2)) cone is the the generating function for the total 
number of branched polymers in the cone. This statement is true regardless of whether 
we put (j) in the bulk and uj near the apex or uj near the apex and (j) in the bulk. On the 



A cone can be thought of as a wedge in the plane, with opposite sides identified. By apex 
angle we mean the angle in the plane between the two sides of the wedge. 



other hand, we see from equation (10) that the exponent ${a) depends on whether or a; 
is near the apex. For if we put uj in the bulk, that is, if we integrate over (2: holding (i 
fixed and near the apex of the cone, we find 

e{a) = 2-— (11) 
a 

while if integrate over holding ^2 fixed we find something different, namely 

e{a) = 2- — -2v (12) 
a 

for a < 2tt. Since we obtain different answers for 9{a) depending on whether we put 
(f) ov UJ near the apex, equation (10) cannot be correct: the correlation function {(f)uj) in 
the cone is not equal to the conformal transformation of (0a;) in the plane, under the 
assumption that (f) transforms as a primary operator. There are in fact two additional 
problems with equations (11) and (12): (1) the the coefficient of ^ is negative, implying 
that the smaller the apex angle, the greater the number of rooted polymer configurations, 
which is impossible, and (2) the exponent 9{a) in equation (12) jumps discontinuously 
from one at a = 27r to about ^ for a just smaller than 27r. Such a jump seems unlikely, 
and in any case is in the wrong direction {9 {a) should increase as a decreases). 

A final point must be considered, namely that the cone is not conformally related to 
the plane, but to the punctured plane. In most theories, the distinction between the plane 
and the punctured plane is irrelevant. A puncture in the continuum theory corresponds 
to a finite size hole in the lattice theory. The presence of the hole changes the interactions 
between sites variables, and so corresponds to the presence of an energy like operator. 
The energy operator on the lattice can be written as a sum of operators in the continuum 
theory. In general, one expects every operator in the continuum theory to appear in this 
sum unless it is forbidden to do so by some symmetry. The lowest dimension operator in 
the sum is the most relevant, and dominates the large distance physics. In most theories the 
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lowest dimension operator with the same symmetry as the energy operator is the identity 
operator, so a defect in the lattice does not change the large distance behavior of correlation 
functions. But in the theory we are considering (and in the Yang-Lee theory) the operator 
(j) has the same symmetry as the energy operator and has a smaller scaling dimension than 
the identity operator. If, therefore, appears in the sum (and there is no symmetry which 
forbids it from doing so) the large distance behavior of correlation functions on a lattice 
with a single defect would be different from the behavior of correlation functions on the 
perfect lattice. In the continuum this means that correlation functions in the plane would 
not be the same as those in the punctured plane. The exponent 9 in the punctured plane, 
or on a lattice with a finite size hole, would therefore be different from 9 in the plane. 
Since correlation functions in the cone are conformally related to those in the punctured 
plane, they would also be modified, as would predictions for 9{a). 

This proposal, while raising interesting questions, has a number of problems. First, it is 
difficult to see how the presence of a finite size hole on the lattice could alter the exponent 
9. To check this, we have enumerated all branched polymer configurations (rooted near the 
origin) with twelve and fewer bonds on a square lattice with the site at the origin removed. 
Branched polymer configurations containing the origin were disallowed. The ratio of the 
number of polymer configurations in the presence of the puncture to the number in the 
plane appears to converge to a constant near .4 as the number of bonds in the polymer 
increases. This indicates that that the defect does not change 9. Second, we expect that 
the presence of the puncture should not affect the boundary conditions at infinity. Thus 
the state at infinity should still be the SL{2, C) invariant vacuum. In the Hilbert space 
formulation, the state corresponding to the operator is an energy eigenstate, and is 
orthogonal to the SL{2, C) invariant vacuum. Hence the partition function Z should not 
be affected by the puncture. Correlation functions, however, will now have an extra 0(0) 
inserted in them. Since we are interested in two point functions, the correlation functions 
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with the extra insertion of 0(0) are three point functions and are computable. Performing 
these computations, one can show exphcitly that the correlation function (0a;) in the 
punctured plane still depends on whether or a; is put at the origin. 

IV. Free Parisi-Sourlas Theory Not Conformally Invariant 

In order to gain insight into what might go wrong with conformal invariance in the 
branched polymer theory, we consider here the simplest theory with a Parisi-Sourlas su- 
persymmetry. This is obtained by setting the potential V{(f)) in equation (1) equal to 
zero. The commuting fields and uj then decouple from the anti-commuting fields i/j and 
ip. Moreover, the fermion part of the free Hamiltonian Hp has the form of an ordinary 
Gaussian model, and is conformally invariant with and t/S transforming as primary fields 
with conformal dimension zero. The remaining part of the free theory is 

Hf = j d!^x[uj{-V'^(j)) + ^uj\ (13) 

or, integrating out the field a;, 

Hp^ j (Pr{V^(pf (14) 

a higher derivative Gaussian theory. For Hp to be scale invariant, must have scaling 
dimension —1, or conformal dimension —1/2. Let us assume, as we did for the branched 
polymers, that (f) transforms under conformal transformations as a primary operator. Then, 
under the transformation z ^ z + e{z), 

ct>{z,z) - (1 - ^)-^(l - '^rh{z,z). (15) 

Taking into account the transformation of the measure (Pr and the derivatives 5^, we find 
the variation in Hp 

SHp = J d\[^^{d(f>){ddcf>) + c.c.]. (16) 

If e{z) is constant, corresponding to a translation, or linear in z, corresponding to a dilata- 
tion or rotation, SHp vanishes identically, while if e{z) is quadratic in z, the integrand is 
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a total derivative, and so 5Hp vanishes in this case as well. These three transformations 
generate the special conformal group, SL{2, C). On the other hand, the integrand is not 
a total derivative for more general conformal transformations. Thus, assuming that (f) is 
primary, the action in equation (14) is invariant under the group SL{2,C), transforma- 
tions mapping the plane to itself, but not invariant under the full conformal group in two 
dimensions. Furthermore, there is no way to transform (f) so that 5Hp = 0. To see this, 
suppose that under 2; — > 2; + e{z) 

with 54> linear in and e. Then 

5Hp = j d^r[(^^^{d(t)){dd(t)) + 2{dd(t)){ddd(t))) + c.c]. (18) 

In order for 5Hp to vanish, we need to choose 5c() in such a way that the integrand is a 
total derivative. Considering the first term in equation 18 we see that this is only possible 
if 5<j) is of the form ae{z)d(l) + hde{z)(l). But we cannot add a term of this kind to the 
transformation law for 0, because this has already been fixed by translational and scale 
invariance (that is to say, if we add term of this form to transformation law for 0, Hp will 
not be invariant under the special conformal group). Thus there is no way to transform so 
that 5Hp = under general conformal transformations. Hence the simplest theory with a 
Parisi-Sourlas supersymmetry is not conformally invariant even at the classical level. This 
may account for the problems encountered in the more complicated theory that describes 
branched polymers. 

V. Numerical Results 

If randomly branched polymers were conformally invariant, the form of the transfor- 
mation law eq.7 for primary operators leads one to expect a a linear relationship between 
the exponent ${a) and the reciprocal of the cone angle a for polymers confined to a cone. 
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A similar relationship would be expected for randomly branched polymers confined to a 
wedge but was not found in a previous analysis of exact enumeration data-^^-^^^. 

We have analysed exact enumeration data for lattice trees confined to a cone. The cone 
is formed by applying cyclic boundary conditions to a wedge, cut from a square lattice, 
in such a way that corresponding lattice sites on the two boundaries of the wedge become 
identified as a single site on the surface of the cone. The number of (weak) embeddings of 
trees rooted at the apex of the cone were enumerated for trees with up to 16 vertices. 

We assume that the number of trees on the surface of the cone diverges as 

Tjv(a) - A^7V-^(")+l (19) 

A is the growth constant for lattice trees on the square lattice and is the same for trees 
confined to a wedge or cone as it is for trees in the bulk lattice'^"*'^'^''*. The yth moment 
of the generating function for tpf{a) will therefore have critical behaviour described by 

G'(y) (x) = J2 Tnx^NV ~ {xc - xf+y-*^ (20) 
N 

From previous work'^^'^*''^^, the value of Xc is known to be 

.,= l/A = 0.19446+''-""l. (21) 
^ ' -0.00002 ^ ' 

To obtain estimates of d{a) a Baker- Hunter confluent singularity analysis^^ was applied 
to the exact enumeration data at each of the cone angles considered. This type of analysis 
was used since we expect the presence of confluences in the generating function for lattice 
trees"^^-^-^^'^. The initial analysis used the central estimate of Xc given above and, since 
Q{a) is bigger than unity for the smaller angles, the first, second and third moments of the 
generating function were analysed^-^^^^. The results are shown in table 1. As this type of 
analysis may be sensative to the value of Xc assumed, we also performed an analysis which 
eliminates the need to specify Xc- To do this the ratio of the number of embeddings in the 
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bulk lattice Tjq and the corresponding number of embeddings in a cone with wedge angle 
a, T/v(ct), was formed 

= TN/TN{a) ~ AT^W"^ (22) 
The generating function for these ratios has a critical behaviour given by 

Gr{x) = J2rNX^ ~ (1 - a;)l+^(")-^. (22) 
N 

This has the advantage that the critical point {xc = 1) is known exactly and it is not 
necessary to use higher moments as the exponent is always greater than 1. The resulting 
estimates of 9 (a) are shown in column 4 of table 1. 

In most cases the spread in the central estimates of 9 (a) in the columns of table 1, for 
a given value of a, is comparable with the error estimate in a single column obtained by 
considering the variation in the estimate from different Pade approximants to the Baker- 
Hunter auxiliary function. In order to make a comparison with the expected linear relation 
we take the overall estimate of 9{a) to have a central value given by an average, over the 
columns of table 1, of the central estimates and error bounds such that all of the central 
estimates fall within these bounds. These values are plotted against 1/a in figure 1. Even 
allowing for some degree of subjectivity in the value of the error bounds the estimates of 
9 (a) are clearly not consistent with a linear relationship. Consequently, although we can 
not altogether rule at short series effects as the reason for this deviation from a linear 
relation, the results reported here together with the earlier work on trees confined to a 
wedge lend support to the conclusion that randomly branched polymers are not conformally 
invariant. 

VI. Conclusions 

In this paper we have presented three arguments against the conformal invariance, at 
criticality, of the supersymmetric field theory that describes randomly branched polymers. 
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We showed first that the most natural formulation of conformal invariance leads to incon- 
sistencies, then that the free Parisi-Sourlas Hamiltonian is not even classically conformally 
invariant, and finally, from numerical ennumerations, that 9{a) is not linear in the inverse 
cone angle a. None of these arguments constitutes a proof; however taken together they 
strongly suggest that randomly branched polymers are not conformally invariant. This is 
interesting, because every other geometrical critical system with which we are familiar is 
described by a conformal field theory. The reason for problems with conformal invariance 
may be related to non-unitarity of the model. Since the theory is non-unitary, the Hamil- 
tonian in the Hilbert space formulation of the problem is not Hermitian, and therefore 
need not be diagonalizable. If H is not diagonalizable, the usual arrangement of operators 
in conformal field theories in highest weight represent aions of the Virasoro algebra^^^"^^*"^, 
would be impossible. 
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TABLE 1 Estimates of 9{a) for varfous cone angles a. Columns labelled 1, 2 and 
3 tabulate results obtained from the first, second and third moment of the generating 
function respectively. The column labelled 4 tabulates the results obtained by taking the 
ratio of the number of trees in the bulk lattice to the number on the cone. 



Angle 


1 


2 


3 


4 


90° 


2.091 


2.134 


2.126 


2.114 




±0.017 


±0.025 


±0.017 


±0.014 


127° 


1.98 


1.96 


1.85 


1.852 




±0.11 


±0.09 


±0.03 


±0.007 


143° 


1.85 


1.83 


1.78 


1.73 




±0.06 


±0.05 


±0.10 


±0.09 


180° 


1.55 


1.551 


1.550 


1.553 




±0.02 


±0.005 


+0.005 
-0.014 


±0.010 


233° 


1.418 


1.4080 


1.39 


1.39 




±0.010 


±0.0007 


±0.09 


±0.04 


270° 


1.2539 


1.250 


1.252 


1.236 




±0.0012 


±0.006 


±0.004 


±0.014 
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FIGURE 1 Overall estimates of 9{a) plotted against 1/a. The bulk value of ^ = 1 
has been included as the value for a = 360°. 
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